Singlet-triplet transition in a lateral quantum dot 



cn 
o 
o 

< 

oo 
(N 



M. Pustilnik,! L.I. Glazman,^ and W. Hofstetter^ 

'School of Physics, Georgia Institute of Technology, Atlanta, GA 30332 
William I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455 
^Lyman Laboratory, Harvard University, Cambridge, MA 02138 

We study transport through a lateral quantum dot in the vicinity of the singlet-triplet transition 
in its ground state. This transition, being sharp in an isolated dot, is broadened to a crossover 
by the exchange interaction of the dot electrons with the conduction electrons in the leads. For a 
generic set of system's parameters, the linear conductance has a maximum in the crossover region. 
At zero temperature and magnetic field, the maximum is the strongest. It becomes less pronounced 
at finite Zeeman splitting, which leads to an increase of the background conductance and a decrease 
of the conductance in the majcimum. 

PACS numbers: 73.23.Hk, 73.63.Kv 72.15.Qm, 
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The Kondo effect in transport through quantum dots 
manifests itself in a dramatic increase of the linear con- 
ductance at temperatures below a certain characteris- 
tic scale T/f (the Kondo temperature). In the simplest 
case Q , a quantum dot behaves essentially as a magnetic 
impurity with spin 1/2 coupled via exchange interaction 
to two conducting leads However, the energy scale 
for intradot excitations is much smaller than the corre- 
sponding scale for real magnetic impurities. Moreover, 
the tunability of this scale in quantum dot devices al- 
lows one to explore various flavors of the Kondo effect 
inaccessible with usual magnetic impurities 

A transition between singlet and triplet states in an 
almost isolated dot was demonstrated Q in a "verti- 
cal" device. At stronger dot-lead tunneling, the conduc- 
tance across the dot has a pronounced maximum at the 
singlet-triplet transition . The maximum can be inter- 
preted 0,0,0 as the Kondo effect with the Kondo tem- 
perature enhanced in the vicinity of the transition. In 
this paper, we investigate the conductance in a "lateral" 
quantum dot device in the vicinity of a singlet-triplet 
transition. The presented below theory of the Kondo ef- 
fect in such devices may help in the interpretation of the 
recent experiments 0, |^ . 

The main difference between vertical dots and lateral 
ones is that in the latter case the number of electronic 
modes coupled to a quantum dot is well defined. A lat- 
eral quantum dot is formed by gate depletion of a two- 
dimensional electron gas at the interface between two 
semiconductors. In this geometry the dot-leads junctions 
act as electronic waveguides. Potentials on the gates con- 
trol the waveguide width, and, therefore, the number of 
propagating electronic modes the waveguides support: by 
making the waveguide narrower one pinches the propa- 
gating modes off. When the very last propagating mode 
nears its pinch-off, the system enters the Coulomb block- 
ade regime. Accordingly, in this regime each dot-lead 
junction supports a single electronic mode [loj . 

Typically, the charging energy of the dot Eq is large 
compared to the mean single-particle level spacing in it. 



5E, which in turn guarantees 2] that 5E is large com- 
pared to Tk- Tk <^ 5E <^ Ec- This separation of the 
energy scales allows us to simplify the problem further 
by assuming that conductances of the dot-lead junctions 
are small. The simplification will not affect the prop- 
erties of the system in the Kondo regime as long as Tk 
remains the smallest energy scale in the problem. On the 
other hand, the coupling between the dot and the leads 
can now be described within the tunneling Hamiltonian 
framework. The microscopic Hamiltonian of the system 
can then be written as a sum of three distinct terms. 



H = H\ 



leads 



^dot ~t~ -^^tunnclingi 



(1) 



which describe, respectively, free electrons in the leads, 
isolated quantum dot, and the tunneling between the dot 
and the leads. With only one electronic mode for each 
dot-lead junction taken into account, the first and the 
third terms in the r.h.s. of Eq. become 

^^loads - Xl^fe^afe'^^afc^' Oi^R,L (2) 

aks 

^tunneling = ^ ] ^cxn'^aks^ns + H.C. (3) 
anks 

A generic model of an isolated quantum dot [the second 
term in the r.h.s. of Eq. |Q] can be written as (l^ 

Hdot = ^ndLdns + Ec(n-N)' - EsS' - BS', 

ns 

where N = J^ns dL^ns and 5 = i XInsV dUa;,,d^^, are 
operators of the total number of electrons on the dot, 
and of the dot's spin, respectively. The parameter J\f in 
Eq. is proportional to the potential on the capaci- 
tively coupled gate electrode and controls the number of 
electrons, N = (N), on the dot. We will assume that 
Af is tuned to a Coulomb blockade valley with an even 
number of electrons N. The third term in Eq. Q de- 
scribes exchange interaction within the dot. Finally, the 
last term in Eq. Q represents the Zeeman effect of an 
external magnetic field, with B being the Zeeman energy. 
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We consider a dot in a state far from the ferromagnetic 
instability the exchange energy Es is small com- 

pared to the mean level spacing SE. Under this condition 
the ground state of an isolated dot with even N is almost 
always a singlet. The only exception occurs when the 
level spacing e = e+i — e_i between the highest occupied 
(n = —1) and the lowest empty {n — +1) single-particle 
energy levels is anomalously small, e ~ Es ^ SE. The 
gain in the exchange energy associated with a formation 
of the triplet state may then be sufficient to overcome the 
loss in the kinetic energy (cf. the Hund's rule in atomic 
physics). The triplet state is formed through a redistri- 
bution of two electrons between the levels n = ±1. Since 
for an isolated dot the occupation S^^d^^ on each 
single-particle energy level is a constant of motion, the 
redistribution must involve tunneling between the levels 
n = ±1 and the leads. At low energies (i3,T <C 5E) 
tunneling to all other (n ^ ±1) energy levels can be 
neglected in the vicinity of the singlet-triplet transition 
(e ~ Eg). Accordingly, in this regime the Hamiltonian 
of the dot Eq. I^J can then be further truncated to that 
of a two-level system with Af — 2 and 



ne/2, 



±1. 



(5) 



In quantum dot systems based on GaAs the value 
of £ can be controlled by a magnetic field H± applied 
perpendicular to the plane of the dot |^ Because 
of the smallness of the electron effective mass, even a 
weak field H± has a very strong orbital effect. At the 
same time, smallness of g-factor in GaAs ensures that 
the corresponding Zeeman splitting remains small 
By linearizing e{H±) in the vicinity of the transition, 
one can make a direct comparison of the experimental 
data with the (calculated) linear conductance across the 
system G = {2e'^ /h)g{e). In addition, one can apply a 
strong in-plane field iJ|| , which would allow the study of 
the conductance dependence on Zeeman energy B. The 
temperature (T) and field (B) dependences of the linear 
conductance, and the bias dependence of the differen- 
tial conductance are qualitatively similar to each other. 
However, the dependence g{e, _B) at T = is easier to 
calculate, and it is this function we address here. 

Tunneling, see Eq. Q, couples the two- level system, 
Eqs. Q and to the two leads. The four tunneling 
amplitudes ta^±i form a 2 x 2 matrix 



In the special case when one of the eigenvalues of t is 
zero, while another one is finite the dot effectively 
interacts only with a single species (a single "channel") 
of conduction electrons. A single channel can screen only 
half of the dot's spin when it is in the triplet state 
Accordingly, the system should exhibit a quantum phase 
transition 'l4'|: the ground state changes its symmetry 



from a singlet to a doublet as e decreases below a cer- 
tain value, £ n ~ Es- The conductance is then strongly 
e-dependent 5(e, 0) cx 9{ec — £)■ At e < ec* (when 
the dot is in the triplet state) the conductance is a mono- 
tonically decreasing function of B, while at e > ec" the 
conductance first increases, and then drops with the in- 
crease oiB 1^. 

In the general situation^ however, both eigenvalues of 
t are finite. Therefore, the dot is coupled to two elec- 
tronic channels, which is sufhcicnt to fully screen the it's 
spin . As the result, the ground state of the system 
is a singlet at all e In other words, when the 

dot is coupled to the leads, the singlet-triplet transition 
turns to a crossover. In order to study this crossover, 
we focus on a special subset 0, 0| of the matrices i 
parametrized as 



t 



1 

71 



■v+i 

V-1 



v+i 



(6) 



Obviously, for v±i ^ both eigenvalues of t are finite, 
so the choice Eq. (|SJ| captures the essential physics of the 
system. 

Since the ground state of the system is not degenerate, 
electrons scatter elastically at T = 0. The amplitudes 
of scattering ^a;aa' of an electron with spin s from lead 
a' to lead a form the scattering matrix §s. The 2x2 
unitary matrix Sg can be diagonalized by a rotation in 
the R — L space to the new basis of channels ti = ± 1 , 



C/S,t/^ = diag {e^'^- } , U ^ e'^^'^'e" 



(7) 



Here r* are the Pauli matrices acting in the R — L space 
(t+ = Tx + iTy transforms L to R). In general, the angles 
i!)q and tpo in Eq. ^ depend on the parameters of the 
microscopic Hamiltonian, and, in particular, on the val- 
ues of £ and B. Here comes the key advantage of Eq. 
with this choice of the tunneling amplitudes, the angles 
are parameter-independent constants: do — 7r/4, ipo — 0. 
Indeed, when written in terms of the operators 

Ipnks ^ ^{nCRks + CLks) , n = ±l, 

[corresponding to — tt/A and (^cq = in Eq. the 
Hamiltonian assumes the form 



nks 



nks 



(8) 

with i/dot given by Eqs. ((SJl. For each n and s this 
Hamiltonian commutes with the operator 



of the total number of electrons with spin s on the "or- 
bital" n. Since N„s commute with each other, the (non- 
degenerate) ground state of H is also an eigenstate of 
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Nns- This implies that single-particle correlation func- 
tions are diagonal in n and s, e.g., {'ipnksW'^li'k' s'i^)) 
Snn'Sss'- Therefore the scattering matrix, which can be 
expressed via retarded single-particle correlation func- 
tions, is diagonal in n and s as well. 

The dimensionless conductance g at T = is related 
to the off-diagonal elements of the scattering matrix by 
the Landauer formula 5=5 X^s l^s ,rl|- With the help 
of Eq. Q the conductance can be expressed via the scat- 
tering phase shifts Sns at the Fermi energy: 



1 
2 



= - sm 



-l,s- 



(9) 



(here we took into account that i?o = 0). The scattering 
phase shifts 5ns in Eqs. lO and Q are obviously defined 
mod TT (that is, Sns and Sns + tt are equivalent). This 
ambuguity can be removed by setting the values of the 
phase shifts corresponding to J\f — —00 (when the dot is 
empty) to 0. With this convention, the phase shifts at a 
finite Af are expressed via the Friedel sum rule in terms 
of the ground state occupation numbers 



ns i^lis^ns)' 



(10) 



Equation IjlUI) is exact in the limit of infinite bandwidth. 
Alternatively, the phase shifts can be extracted directly 
from the finite-size spectra obtained by the numerical 
renormalization group [l^ . 

In the Coulomb blockade valley the number of electrons 
on the dot is fixed: N — J2ns ^ns = 2. Accordingly, the 
phase shifts satisfy 



2tt. 



(11) 



Additional relation for the phase shifts follows from the 
invariance of the Hamiltonian (^-Q with respect to the 
transformation {s,B) —^ (—s,—B): 



Sns{B) = Sn,-s{-B). 



(12) 



Consider first the limit of zero field, B = Q. From 
Eqs. (dJ and it follows that 



(13) 



At the triplet side of the crossover e <^ Es both levels 
in the dot are singly occupied, so that Eq. ((Tn|l yields 
Sn = 7r/2. On the contrary, at e ^ Es the level n = —1 
is doubly occupied, while the level n = -1-1 is empty, 
which corresponds to (5_i = tt, S+i = 0. When e is 
tuned through the crossover, the difference of the phase 
shifts. As = (5+1 — (5_i [see Eq. monotonically de- 
creases from to — TT. It follows then from Eq. that 
the dependence of the dimensionless zero-field conduc- 
tance (7(e, 0) on e is nonmonotonic. The conductance 



reaches its maximum g = 1 at some e = Eq correspond- 
ing to As(eo) = — 7r/2, and falls off monotonically with 
the distance |e — eo| to this point, see Fig. The en- 
ergy £0 — 2£'5 may be identified with the center of the 
crossover region. The conductance 5 ~ 1 within this re- 
gion. Later on we will relate the width of the crossover 
region Tq to the parameters of the Hamiltonian 

The effect of the Zeeman energy B on the conductance 
accross the dot depends on how far its parameters are 
from the crossover point. In order to study the influence 
of a small B, one can expand the phase shifts in a se- 
ries 0. Taking into account Eqs. (QU, (dJ, and ifT^ . 
we obtain 

Snsie, B) = Sn + s {B/Tn) + u [B/T'f + 0{B^), (14) 

where n = ±1, s = ±1, and Sn, Tn, and T' depend on e. 
The sign of the linear in B term in Eq. (|14|l is fixed by the 
following argument. The spin polarization Nn^+i — Nn.~i 
of electrons on each level (n = ±1) grows proportionally 
to B. It follows then from Eqs. (^01 and that r„ > 
for n = ±1. In the next order in B, magnetic field favors 
a triplet over the singlet state of the dot. Therefore, the 
difference iV+i,s — ^-i,s increases with \B\, which fixes 
the sign of the second order in B contribution in Eq. H14|l. 




FIG. 1: Dimensionless conductance across the dot at T = 
and B — 0. The conductance decreases with _B at e_ < e < 
e+ and increases with B outside this interval. 

Substitution of Eq. lfTI|l into Eq. (O yields the low-B 
asymptotics of the conductance, 

g{e, B)^g + (1- 2g) {B/Bof + 4^5(1-5) (B/T'f , 

(15) 

with g ~ g{e, 0) and 



(16) 



So(£) 



En 

n— ±1 ^ ' 



Around the crossover point g ~ 1, the last term in 
Eq. H15|) vanishes, and 

g{e,B)^l-{B/B^f , |£-£o|«ro. (17) 

Away from the crossover (5 ^ 1), Eq. I|15|l yields 

g{e,B)^g{e,t)) + {B/Bof , |e-eo|»To. (18) 
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It is clear from Eqs. H17() and (|18|l that the conduc- 
tance varies with the magnetic field in the opposite di- 
rections in the vicinity and far away from the crossover 
point. It should be emphasized that this is a generic 
property rather than a consequence of the specific choice 
of the model, Eq. ©. However, the precise borders of 
the crossover region e — e±, see Fig. as well as zero- 
field conductances g{e±,0) at these points are model de- 
pendent. Note also that Eq. (|15|l is applicable only at 
B <C Bo{e). The conductance at higher fields is a non- 
monotonic function of B, see Refs. UlIBIlllia. 

The Hamiltonian Eq. (jSJ is identical to that employed 
previously 0] to study transport through a vertical dot. 
Accordingly, the thermodynamic properties of a lateral 
dot coincide with those of a vertical dot. However, 
the existing experiments test transport properties rather 
than thermodynamics. The electron current operators 
are different for the two models, hence the dependences 
g{e,B) are different as well. For instance, the conduc- 
tance through a vertical dot at T = is large (~ 4e^/ft,) 
at the triplet side of the crossover and decreases with B, 
in contrast with Eq. (|18() . 

At |e — eol ^ To the parameters g{£,0) and Bo{e) 
entering Eq. (|18|l can be estimated with the help of the 
perturbative renormalization group 0- In this regime 
the e-dependence of all observable quantities is governed 
by the parameter 

x{e) - (e - eo)/To, 1 « |a;(£)| « SE/To. (19) 

The estimated width of the crossover region satisfies 

H5E/To)^ 7 = k+i/«-i|, (20) 

where v is density of states in the leads. The function 
(7(7) has a minimum C ~ 0.36 at 7 = 1 and goes over 
to 1 as 7 ^ 0, 00. The scale Tq also plays the part of the 
Kondo temperature as it determines the T-dependence 
of the conductance in the crossover region: the peak in 
g{e) disappears at T ^ Tq- Typically, Tq is larger than 
the Kondo temperatures in the nearby Coulomb blockade 
valleys with odd number of electrons on the dot 6, 7]. 

The zero-B conductance entering Eq. H15|l away from 
the crossover point (at \x\ ^ 1) is 



3(e,0)cx(ln|x|)- 



(21) 



Finally, the characteristic magnetic field Bo{e) defined 
in Eq. (|16|) has different asymptotes at the triplet and 
singlet sides of the crossover: 

r |2;|-^(ln|a;|)2V(A+i)^ 2; < 0, 
Bo/To « <^ (22) 
[ :E(lna;)2(^+2)/(^+i), x>Q 

with A = 2 -I- \/5 w 4.2 and \x\ > 1. 



The dependence of the differential conductance dl/dV 
on the source-drain bias V at eV ^ T,B is qualitatively 
similar to the dependence of the linear conductance on 
B at B ^ T. The similarity stems from the fact that 
dl/dV is determined by the electron transmission coef- 
ficient at energies w ~ eV, where co is measured from 
the Fermi level 0. The dependences of the transmis- 
sion coefficient on to and B are controlled by a single en- 
ergy scale: the a; = limit of the transmission coefficient 
governs [2j 'is'l the linear conductance g{e,B). There- 
fore, Eq. (|17|) implies also that dl/dV decreases with V 
in the crossover region shown in Fig. ^ On the con- 
trary, away from the crossover the conductance is small, 
dl/dV - {2e^/h)g{e, 0), in the domain F < - B^/e 
and is increasing with V , cf. Eq. (|18|l . The "window" 
is narrow at the triplet side of the crossover but broadens 
approximately linearly with the distance to the crossover 
at the singlet side of it, see Eg. 1I22II . This is in agreement 
with numerical simulations |l8l | . 

The zero-bias suppression of dl/dV was observed in 
recent experiments |9|] . A measurement in the limit of 
very weak tunneling, necessary for the proper character- 
ization of the spin state of the dot, was not possible for 
the studied devices. However, the observed |3 asymmet- 
ric behavior of AV across the point where the conduc- 
tance has a maximum is precisely the expected behavior 
in the vicinity of the singlet-triplet crossover. The asym- 
metric behavior of AV^ across the crossover may explain 
the observed 9] splitting of the Kondo peak in a certain 
range of gate voltages. 
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